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We study the quantum dynamics of an optomechanical setup comprising two optical modes and 
one mechanical mode. We show that the same system can undergo a Dicke-Hepp-Lieb superra- 
diant type phase transition. We found that the coupling between the momentum quadratures of 
the two optical fields give rise to a new critical point. We show that selective energy exchange 
between any two modes is possible by coherent control of the coupling parameters. In addition we 
. . . also demonstrate the occurrence of Normal Mode Splitting (NMS) in the mechanical displacement 

spectrum. 

O ■ PACS numbers: 42.50.-p,42.50.Ct,42.50.Wk 

(N ■ 

(D ■ I. INTRODUCTION 

\^ ' Electromagnetic coupling between optical cavities and nano-mechanical resonators has led to new quantum mechan- 
: ical behaviour of macroscopic quantum systems [M. Substantial experimental progress has been made to realize such 
(-H novel quantum systems at the single photon level [2ru] • At the same time several theoretical works have appeared in 
O ^ the literature that study a single optical mode coupled to a single mechanical mode, in the strong coupling regime. 
^ I , Quantum effects are found in these hybrid systems if the coupling between the light field and the mechanical resonator 
(~| becomes much larger than the cavity decay rate and the mechanical oscillation frequency. Quantum nonlinearities 

can also be introduced into the system by placing a Kerr medium inside the optical cavity jS] . 
^ ^' It has been demonstrated experimentally that by coupling two optical modes to a mechanical oscillator, the quantum 
nonlinearity can be enhanced significantly jol-fl^. This enhanced nonlincarity has potential application in quantum 
nondemolition (QND) measurement of phonon number p^ . Besides QND measurements, optomechanical systems 



have potential applications in quantum information processing jl4| . The kerr type nonlinearity can form the basis of 
K*" ' all optical switch with application to engineer a quantum phase gate for photonic or phononic qubits. The mechanical 
' mode can also serve as a quantum memory 15]. The interaction between the optical and mechanical degrees of 
• freedom gives rise to a quantum interface between solid-state, optical and atomic qubits [16i] . Optical modes in 
quantum information processing units can transfer information over long distances and on the other hand acoustic 
excitations (phonons) can store information for an extended period. A hybrid architecture composed of optical modes 
■ and mechanical modes can thus be fruitfully utilized to design a quantum communication and quantum information 
, processing unit which can store and transfer information coherently. Theoretical studies related to two mode system 
have been performed much in detail (T7l-l20| . 

Given these new developments in this field, in this work, we show that the dynamics of the linearized Hamiltonian of 
the two optical modes coupled to one mechanical mode is equivalent to that of the Dicke model [2l| . We analytically 
calculate expressions for the normal modes of the system and show that the system can undergo a Dicke-Hepp-Lieb 
^ . type superradiant phase transition by varying the various coupling constants. We demonstrate the possibility of 
selective energy exchange between any two modes. We also study the squeezing variances of the three modes near 
the quantum critical point. We demonstrate the normal mode splitting (NMS) in the displacement spectrum of the 
mechanical oscillator. 

II. DICKE-HEPP-LIEB SUPERRADIANT PHASE TRANSITION MODEL 

The optomechanical setup that we consider here (Fig.l) is composed of an optical cavity with two optical modes (de- 
noted by operators ai and 02) and one mechanical mode (denoted by 6). This system is described by the Hamiltonian 

mm 

H = Ho + Hint, (1) 



where 



Mechanical mode 



Figure 1: Schematic representation of the hybrid optomechanical system consisting of two optical modes denoted by ai, 02 and 
one mechanical mode b. Also shown are the coupling rates between the various modes. 



Hq = hLL!ia\ai + huJ2a\a2 + hajmb^b, 



(2) 



Hi, 



hgia\ai + hg2a\a2 — hG{aia2 + a2o\) {b + b^) 



(3) 



Here Hq is the bare energies of the two optical modes with frequencies Ldi,2 and that of the mechanical mode with 
frequency ujm- The term Hint denotes the interaction between each of the two optical modes and the mechanical 
oscillator with coupling rates 31,2- The coupling between the two optical modes via the mechanical oscillator is G. 
This kind of Hamiltonian is found typically in a setup with a membrane in the middle of the cavity 1^ . The interaction 
Hamiltonian results from the space dependence of the cavity mode frequencies. 

We are interested in studying the dynamics of fluctuations of the system. To this end, we rewrite the operators 
ai^2 and b around their mean values as, ai^2 — ^ cn,2 + 01,2 and 6 — /? + 6. Here ai^2 and /3 are the mean field values 
of the optical modes and the mechanical mode respectively whose values are given in the Appendix I. As a result, we 
get an effective Hamiltonian for the system after retaining terms which are bilinear. 



H = SfJiajai + ?iJ72a2'^2 + fujJmb^b — hGi{a\ + ai) + hG2{a\ 



- 0,2) 



{b + b'l) + h\{aial + a2a\), 



(4) 



where fli.2 = wi_2 — 2j3gi^2, Gi^2 — 5i,2Cki,2 ~ Ga2,i and A — 2Gf3. 

The Hamiltonian Eqn.(4) is bilinear in the fiuctuation operators and can be easily diagonalized. This is accomplished 
by the following position and momentum operators for the three bosonic modes. 



--{al + ai), px i 



Oi) 



(5) 



y = 



+ 02), Py = i\h^{a\ - a2) 



(6) 
(7) 



The Hamiltonian (after ignoring constant terms) in terms of these operators is written as 



nlx^ +pI + nly^ +pI + ujIz^ +pI- 4Gl^/n^xz - 4G2Vih^yz + 2A (^ih^hxy + ;^^=) 
We now rotate the coordinate system in the following way: 



^=2 



• (8) 



X = g'i(cos7i + cos 72) + 92 sin 71 + 93 sin 72, 



(9) 



y = -qi sin 71 + 52(00371 + cos 73) + qs sin 73, 



(10) 



z = -qi sin 72 - 92 sin 73 + 93 (cos 72 + cos 73) 
where the angles 71, 72 and 73 are given by 



(11) 



tan 271 = —2 — -2-, tan 272 = —p^ — -1^, tan 273 - 



(12) 



This rotation eliminates the xy , xz and yz interaction terms in the Hamiltonian. In order to eliminate the PxPy 
interaction term, we further make the following transformation, 

Px = Pi cos P +P2 sin P, Py = -pi sin /S +P2 cos /?, (13) 
where the angle /3 is given as cos 2/3 = 0. The Hamiltonian then takes the form of three uncoupled oscillators. 



_ 1 J 4q! , , elql o 2,2] 

^ ~ 2] ~~2~ + P3 I" > 



where 



(14) 



(15) 



(16) 



(17) 



epi = 1 



A 



=5 — 1 



A 



We now requantize H by introducing the following three bosonic modes. 



(18) 
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93 = 



C3 



P3 = M/ TT 



C3 



(21) 



We arrive at the final diagonal form as, 



H = exclci + eY4c2 + ezcjcs + -^[(^x + + ez] (22) 
where, ex = Cxi/e^, er = Cy ^^"^ = ^z- 

It can be noted from equations [T51 [TOl [T71 and that the value of the eigenvalue ex corresponding to the c\ci 
bosonic mode can become imaginary when e^ep^ < 0. This happens when A > v^TTlT^- This means that the system 
undergoes different bahaviours depending on being in the "normal phase" i.e when e^Cpj > or in the so called 
"super-radiant phase" when e^Cp^ < 0. The eigenvalue ey corresponding to the c\c2 bosonic mode remains always 
real. The eigenvalue ez can also become imaginary when we scan the coupling constants Gi or G2. Since we are 
interested only in the normal phase regime, we will not go into the details of the phase transition regime. The model 
described above is valid only in the normal phase. The excitation energies ex (thin solid line), ey (dashed line) and 
ez (thick solid line) of the 3-mode Dicke type Hamiltonian of equation!^ as a function of coupling parameter A (plot 
(a)), Gi (plot (b)) and G2 (plot (c)) is shown in Fig.l. The ex.y are the optical branches while ez is the phononic 
branch. Selective energy transfer is seen in these plots. In plot (a), energy transfer between the ci optical mode and 
the C2 optical mode is observed while the energy of the C3 phononic mode does not change. The eigenvalue ex displays 
the phase transition at the critical coupling constant Ac. The critical coupling constant is Ac = •\/^i^2. The critical 
point at Ac is a consequence of the coupling between the momentum quadratures of the two optical modes. This is a 
new feature that we observed in this system. Interestingly we found that for some parameters the eigenvalue ex takes 
negative values, signifying instability. This happens when both gp^ < and < 0. The value of Xus above which the 
system is unstable is: 



VnTn^ ■ ^^^^ 

In absence of any momentum coupling, A„s is the usual critical point Ac- For plot (b) energy exchange between 
the ci and C3 modes is seen. Now in this case ez demonstrates the phase transition at a certain critical coupling 
parameter Gic- The value of Gi above which the excitation energy ez becomes imaginary is: 

J {ill + v„ - ^{^i <r + i6Gin2iOr..r - i^i ^ir , , 

Gi > ^ 77n= ■ (2^) 

The superradiant phase transition in the phonon excitation spectrum indicates that complete energy transfer from 
the mechanical mode to the optical mode occurs at the quantum critical point. Similarly in plot (c), energy exchange 
between C2 and C3 modes is clearly visible. 



III. MODE SQUEEZING 

In this section, we investigate the squeezing behaviour of the three bosonic modes. A bosonic mode is said to be 
squeezed if the uncertainty in either of its position or momentum quadrature is less than the uncertainty in a coherent 
state [22|. A coherent state is a minimum uncertainty state which satisfies (Aa)^(Apct)^ = 1/4 [a — x,y,z] and 
with uncertainties equally distributed between the two quadratures. A bosonic field is squeezed if (Aa)^ < 1/2 or 
(Apa)"^ < 1/2 [23. The two quadratures variances of the original field modes ai, 02 and b are defined as (Aa)^ =< 
> — < a >^ and {Ap^)'^ =< Pa > ~ < Pa The relationship between the original bosonic modes [oi, 02, b] and 
the transformed bosonic modes [ci, C2, C3] are given as, 
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Figure 2: (Color Online) The excitation energies ex (thin solid line), ey (dashed line) and ez (thick solid line) of the 3-mode 
Dicke type Hamiltonian of equation [22] as a function of coupling parameter A (plot (a)), Gi (plot (b)) and G2 (plot (c)). The 
ex,Y are the optical branches while ez is the phononic branch. Selective energy transfer is seen in these plots. In plot (a), 
energy transfer between the ci optical mode and the C2 optical mode is observed while the energy of the C3 phononic mode 
does not change. For plot (b) energy exchange between the ci and C3 modes is seen while in plot (c), energy exchange between 
C2 and Ci modes is clearly visible. 



1 ^ cos 71 + COS 72 

2 ^Vtiex^/f^ 



+ 



sm 71 



sin72 



4(1^1 +e^)+C3(0i-e^)j} 



(25) 



6 



4 = ^{ 



- sm 7i 



+ 



+ 



COS 7i + cos 73 



sin73 



(26) 



+ 



sm 72 
sin 73 

COS 72 + COS 73 



One can trivially show that the various variances are, 



A{^in + e^) + C3(Wm - ^z) } 



(27) 



(\ ^2 1 /i , (cos7i +cos72)% ^ — . sin27i ^ — sin^ 72 



2^1 



(28) 



n2 1 \^^^'^^^^( I — n ^ , (cos 71 + cos 73)2 , — , sin2 73 

\^V) = T^FT i 1 H (v^^2 - Cx) H (v^^2 - Ey) H (1^2 - e^) 

I ej; €2 



(29) 



/A \2 1 /1 , '('2/ ; — \ , sin 73/ , — X , (COS72 +cos73)^ 



(Wm - e^)| 



(30) 



(Ap^)2 ^ ^ L ^ (^os7i + cos 72)2 ^ , sin^ 71 



-(ex - -v/e^^i) + 



v^. - + ^^(e^ - ^1) \ (31) 



(Ap,)^ = ^ ^ 1 



sin^ 71 



pi"2 



(fix 



pii^2 



(cos 71 + COS 73)^ 



P2"2 



(ey - V^^2) 



sm 73 

02 



(e^ - f^2) 



(32) 



(Ap,)^ = =^ 1 



sin^72 . , — X. sin2 73 — ^ ^ (cos 72 + 00^73)2 1 



/em W„ 



Fig. 3 displays the squeezing variances corresponding to Fig. 2. In Fig. 3(a), the momentum variances {Ap^)^ and 
{Apy)^ are unsqueezed and diverges (more unsqueezed) as A approaches Ac. On the other hand the position variances 
(Aa;)^ and (Ay)^ are already squeezed and become more squeezed as A approaches Ac. In contrast, the Figures 3(b) 
and 3(c) illustrates that the momentum variances which are initially unsqueezed, appoaches the value 1/2 and becomes 
slightly squeezed as the coupling parameters approach the critical value. The corresponding position variances which 
are initially squeezed diverges and becomes Lmsqucczcd as the critical point is approached. Thus we conclude that the 
quadratures which arc squeezed are best suited for quantum measurements since squeezed quadratures have reduced 
quantum noise i.e, if we make quantum measurements while varying A, then it would be more appropriate to make 
measurements of the position quadratures. On the other hand, while varying Gi or G2,the position quadrature is a 
better candidate as long as it squeezed. Near the critical point, measurements on the momentum quadrature would 
be better suited. 
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G2/C0n 



Figure 3: The squeezing variances corresponding to Fig. 2. (a): Plot of (Aa;)^ (dashed hue) ,(Aj/)'^ (dashed hne), [Ap^]^ (solid 
line) and (Apy)^ (solid line) as a function of \/uJm for Qi = ujrn, ^2 ~ i^m, Gi = O.OlcJm and Q2 = O.OlcJm. Note that (Aa;)^ 
and (At/) are coincident and the same for [Ap^f and [Apyf. (b): Plot of (Aa;)^ (dashed line) ,{Azf (dashed hne), (Ap^ ' 



(solid line) and (Ap^) (solid line) as a function of Gi/ujm for Q,i = lj„ 



X = O.Olw™ and G2 = 0.01a;„. (c): Plot 



of (Ay) (dashed line) ,{Az) (dashed line), (Apy) (solid line) and (Apz) (solid line) as a function of G2/uim for Q.\ = a;„ 
^2 ~ i^m, A = O.Olcjm and Gi — 0.01a;m 



IV. DISPLACEMENT SPECTRUM: NORMAL MODE SPLITTING 



Here we show that coupling between the mechanical mode fluctuations and the two cavity mode fluctuations leads 
to a splitting of the normal mode into three modes (Normal mode splitting (NMS)).The optomechanical NMS however 
involves driving three parametrically coupled nondegenerate modes out of equilibrium. The NMS does not appear in 
the steady state spectra but rather manifests itself in the fluctuation spectra of the mirror displacement. 

Using the Hamiltonian given by Eqn. ([?]), we get the linearized equations of motion which are represented as follows: 
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di{t) = [iAi - ^]ai{t) + iGiibit) + bHt)) ~ iXa2{t) + V7^a,„,(i), 
d2{t) = [iA2 - ^]a2(<) + iG2{b{t) + b^t)) - iXai{t) + ^a,n,it), 



(34) 
(35) 



b{t) = -{iLj^ + jrn)b{t) + iGi{ai{t) + a\{t)) + ^62(02(0 + 4W) + V^^(t), (36) 

where Ai — —Qi, Gi ~ giai — Ga2, A2 = —Q2 and G2 — 92012 ~ Gai. Now, we rewrite the above equations 
in terms of amplitude and phase quadratures for the system with Xi{t) = [ai{t) + al(i)], Yi{t) — i[a\{t) — ai(i)], 
X2(t) = [a2{t) + 4{t)], Y2{t) = i[aUt)-a2{t)], Q{t) = [b{t) + bHt)], P{t) = t[bHt) - b(t)], X„,,{t) = [a„, (<) + (i)], 
y„,(i) = i[aj„^(i) - ai„,{t)], X„,(<) = [ain^{t) + a|„^ (i)] and Yin^{t) = i[a|„^(i) - atn^it)] as: 



Xi{t) = -AiY,(t) + \Y2{t) - ^Xi{t) + ^X.nAt), 



Yi{t) = A,X,{t) + GiQ{t) - \X2{t) + ,/l^^n^{t) " ^^(0, 



X2{t) = -A2Y2{t) + AFi(i) - ^X2{t) + Vl^X,^,{t), 



Y2{t) = A2X2(i) + G2Q{t) - XX lit) - ^Y2{t) + V7^r™,(0, 



(37) 
(38) 
(39) 
(40) 



Qit) = iJrnPit), 



(41) 



Here W{t) 



P{t) = -UJ,nQ{t) - J,nP{t) + GiXiit) + G2X2{t) + W {t) . 

/lm{£,Ht) — ^{t)) satisfies the fohowing correlation [23j : 



(42) 



1 + coth 



\2kBTj 



where fcs is the Boltzmann constant and T is the temperature of the reservoir. 

The displacement spectrum for the mechanical mode of the system is obtained by using 



(43) 



where 



Qiio) 



1 

47r 



Ai{uj) + A2{uj) + A^juj) + Ai{uj) + Ar^jio) 



(44) 



(45) 



using the correlation functions given in the Appendix II. The values of Ai{uj), 74.2(0;), 74.3(0;), ^44(0;), ^15(0;) and 
B(uj) are given in the Appendix III. 

Fig. 3 shows the displacement spectrum {Sq{uj)) as a function of dimensionless frequency (o;/wm)- The NMS is 
associated with the mixing between the mechanical mode and the fluctuation of the two cavity field around the steady 
state. Clearly we see three modes in the displacement spectrum indicating a coherent energy exchange between the 
mechanical mode and the two optical modes. An important point to note is that in order to observe the NMS, the 
energy exchange between the three modes should take place on a time scale faster than the decoherence of each mode. 
On the positive detuning side, the observation of NMS is prevented by the onset of parametric instability. 

The experimental prospects for various parameters used in the main paper are illustrated as follows. Mechanical 
frequency of optomechanical system can take the value 2tt x 73.5 MHz with corresponding damping rate 2tt x 1.3 KHz 
[2^ . High Finesse optical cavities can decay with a rate nearly 27r x 10 MHz [2^, [2^. Optomechanical crystal setups 
record a coupling rate of 27r x 1 MHz [l^. Such coupling rates can be enhanced above 10 MHz by using nanoslots 
[2^ which increases the local cavity field in these structures. The limit huj„i << fi'jm « k^T is always taken into 
account for various optomechanical experiments [29143^ . 
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Figure 4: Plot of displacement spectrum {Sq{u})) as a function of dimensionless frequency(aj/ctj„i). Parameter values are : 

7m = W~'^UJm, Al = — l.ScJm, A2 — — l.SoJm, = 0-2uJm, 7c, — 0.6LOm., G = 1.5cJm, /? = 0.06, Gl = 2aJm , G2 = GtOm aud 

kBT/hw^ = 10^ 

V. CONCLUSION 

We have presented how the existence of a quantum phase transition in a hybrid optomechanical system can be 
exploited to selectively transfer energy between two optical modes and a single mechanical mode. By scanning one of 
the three coupling parameters, we can selectively open one channel to transfer energy between any two modes while 
at the same time close the other two channels. At the quantum critical point a complete transfer of energy is possible 
between any two modes. The study on the momentum and position variances tells us that the quadratures which 
are squeezed are best suited for measurements due to less amount of quantum noise. The normal mode splitting 
shows distinct three peaks which again demonstrates that the process of coherent energy exchange between the three 
modes is possible. The fact that phononic modes can store energy for a longer duration and photonic modes can 
transfer energy over long distances makes this hybrid optomechanical system useful in next generation of quantum 
communications and quantum information processing units. 

VI. APPENDIX I 

Starting from the Hamiltonian of Eqn.(IT]), we derive the Heisenberg equations of motion for the operators ai(i), 
a2(t) and b{t). This yields. 



ai 



hit) = -lUJiaiit) + igiai{t){h{t) + b\t)) - iGa2{t){b{t) + h\t)) - ^ai{t) + V7^a„,(i), 



(46) 



d2{t) = -iu)2a2{t)+ig2a2{t){b{t)+b\t))-iGai{t){b{t)+h\t)) - 2^a2{t) + ^a^n^t), (47) 



bit) 



Mt) + i9ia\{t)a^{t) + ig2al{t)a2{t) ~ tG{ai{t)al{t) + a2{t)a\{t)) ~ -i^b{t) + ^^{t). 



(48) 



where Uimit) and ain^it) are the input noise operators for the two optical modes, ^(t) is the noise operator arising 
from the brownian motion of the mechanical mode. Now, we find the steady state parameters for the different 
operators by factorizing the non-linear algebraic eqs. (|46p - (j48p and setting their time derivatives to zero, given as: 



iGa2{P + P^) 



[igiiP + P^) - {lu^ 



-)]' 



(49) 
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0.-2, 



iGai(/3 + /3t) 



[»52(/3 + /3t)-(»L.2 + 2|^)]' 



(50) 



(51) 



(iWm + 7m) 

where ai, 02 and /? are the steady state values for the two optical cavity modes and the mechanical mode respectively. 



VII. APPENDIX II 

The correlation function for the noise operator arising from brownian motion in the fourier space is given by [SS 



{W{uj)W{J)) = 2^w^ |l + coth { ) 1"5(^ + ^')- 



(52) 



The correlation functions for the amplitude and phase quadratures of the various input noise operators in the fourier 
space are given as follows [33j : 



(X,„,(w)y„,(c^')) = (X^„,(c^)r,„,(w')) = 2z^,5(w + w'), 



(54) 



(K,„,(c<.)X,„,(c^')> = (K,„,(c<.)X,„,(c^')) = -2^7r(5(c^ + c^')■ 



(55) 



VIII. APPENDIX III 

The coefficients in equation (j45p are given as follows: 



A^(lS) = W^(tj)w„ ( A2 - ^2 ^ ^ + iw7,., ) Ci(w)C2(w), 



(56) 



A2(c^) = X,M^. ( - ^ ^ + ico7e, ) Cx{^)Cz{u:), 



(57) 



^3(0;) = Fi„i (iS)^f^^ ( A2 - + ^ + za;7c2 j Aw,„G2 (^zw + ^ ) C2(w) 



A2A2 - Ai ( A2 - ^ 2i + ^^^^^ ) J, C'3(c^) 



(58) 



^4(w) = X„^^[uj)^/^_ ( A2 - W 



- ia;7c. 



V„G2 {Ai (^c. + + A2 + I 
+w™G2G2(c^) (*^ + ^) +Aw™GiGi(l^){Ai (zc^ + :^) +A2 



(59) 
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XCsico) jA^Ai - A1A2 (^Ai - + ^ + zc7,,^ + Ciico) {ico + ^) } 

/ 2 



(60) 



where 



(61) 



(62) 



C2(a;) = ( Ai-a;2^^+iu;7e,) (A?-a;2 + ^+iw7e, ) +A* 



+ 2A2{(iu;+^) (iu; + ^)-AiA2}, 



(63) 



+ Au;mG2Ai 



AwmG2A2 ^iw - 



2 y ' 



(64) 



C4(a;) =(^Al-oj'' + ^+ icvy^,^ ^2(0;) [(w^ - + iuj'y^)Ci{uj) + LOmG'' A2 (iw + ^) 



-Aw„GiG2 (zw + ^) 



(65) 



G5(a;)=G3(w) 



AG2Gi(w) 



AG2A2 + Gi ( A^ - ^2 ^ ^ + ia;7,. 



)} 



Ai ( A2 - + ^ + iu;7e. 



A2A2 



(66) 
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